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its logarithm is a normally distri- 


buted variableat eachstress level. Eachstress level is individually ana 

lyzed by assuming the sample moments to be good estimates of population 
moments. This procedure often leads toinconsistent results in the relative 
amounts of scatter at various stress levels Sometimes a high stress is 
shown to have more scatter than a low stres in such analyses The life 
to stress relationship consequently becomes obscured, and most analysts 


work on the assumption thata linear relationexists betweenlogN and logS 
em, in spite of the inconsistent arrangement of 
points brought about by classical statistical analysis. While a linear re- 


lation between log N and log S witha fixed s! ypeisa usable approximation 
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alid at all stress levels and which will tell us how much scatter can be 
expected at any stres Most experts have concluded that the line of best 
performance inan S-N diagram is less steep than the line of worst per- 
lormance Why Because as they say t s only reasonable to have 
more scatter at lower stresses than at higher stresses, and hence the tw 
lines mustconverge as we move upward on the stress axis But, do they 
actually converge toward each other in a linear fashior and come toa 
point of intersection Furthermore is there any relationship between 
the rate of convergence andthe endurance limit? Theseand other questions 
cannot be answered until we firstestablish the fundamental relation between 
stress, life and probability of failurs 
The present work is an attempt to arrive at this relation in an axio- 
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DEFINITIONS AND ASSUMPTIONS 


DEFINITION I (q-life): The q-lifeof a stressed specimen is the number 


of cycles for which the probability of fatilure is q (Denoted by N, } 


DEFINITION II (Characteristic life): The characteristic life ofa stressed 





specimen is the number of cycles for which the probability of failure is 
1 - l/e 0.032 { De noted by 0.) 
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DEFINITION III, (Minimum life): The minimum 


fe ofastressed specimen 





is the minimum number of cycles which any such specimen will endure 


under the given stress Some authors call this the ''Sensitivity limit 


(Denoted by Na, } 
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DEFINITION IV, (Static strength): The static strength of a specimenis 


that stress at which the specimen has a minimum life of zero cycles and 





n a characteristic life of one cycle (Denoted by U) 


S DEFINITION V, (Endurance limit): The endurance limit ofa specimen is 





e that critical value of stress at or below which the specimen will never 
S fail under repeated loading, i.e., the stress under which the specimen 
S has a characteristic life of infinity (\Denoted by L). 


Using the above definitions, we now proceed to build a theory of S-N 
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We shall justify assumption I, as well as arrive ata formula for (gq), in 


™ appendix l. This assumption essentially staté that a life distribution 
is } functionat a given stress 1S determined bDby exactly two parameters, 1.e., 
N and 8 L (q) being the same for a! stress levels 
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ASSUMPTION II: Any S-N equationmust be dimensionally homogeneous. 


This assumptionis commonly adoptedas an axiomin dealing with relations 
1 i I 





} between physical quantities, and we adopt it here as well 
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, THE BASIC VARIABLES OF THE GENERAL S-N EQUATION 

d 

is |} Every fatigue specimen of a given material has a definite endurance 
limit, L, and a static strength, U. The magnitude of the applied stress 
shallbe denoted byS. Therefore, there are three stress variables enter- 

en ing into the S-N equation, i.e., L, U, andS 

re 


The S-N equation is an expression for the q-life at stress 5S, i.e ., a 
formula forN, . interms ofthe stresses L, U, andS, andother necessary 
variables. These other necessary variables are determined from as- 
sumption I, which states that Ng sc is a function of the three variables 


N O(q), and 6 Since the S-N equationaiready containsL, U, andS, 
” _ ¥* , Ss : 





we are permittedto take only two of the three variables No 
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assumption I in that it implies that O(q) is a function of L, U, and S. We 


therefore take either of the other two possible pairs from No, sg, Q(q), and 


8. as our additional variables. Let us take the pair Ny, 


the S-N equation becomes 
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Hence, (1) and (2) imply that 6. is a function of L, U, and S, which isa 
reasonable conclusion. We therefore adopt the form (1) asa permissible 


S-N equation. 
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Therefore, the S-N equation (1) may be written as follows: 


THE CONSEQUENCES OF ASSUMPTION II 
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In terms of length, force, and life, the dimensional matrix of this set of 
variables is 


Ng,s L U S No. s Ng, U 
length 0 -2 +2 -2 0 0 
force 0 ] ] ] 0 0 
life ] 0 0 0 l ] 


The rank ofthis matrixistwo. Therefore, thereare 6 
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A particular solution, satisfying these two conditions, is 


N N ‘NM... = =] , (6) 


which may be taken asa generalS-N equation, provided that itcan be made 
to fit experimental data At the present time, there does not exista suf- 
ficiently wide collection of fatigue life data from independent sources to 
enable us to verify (6) very thoroughly. However, a critical examination 
of (6) indicates that it is very flexible, and can be made to fit very closely 


any fatigue situation which may come up for consideration. 
APPENDIX I 
. ; * 
Justification of assumption I: A.M. Freudenthaland E. J. Gumbel” have 
shown that ata givenstressS the fatigue life of a specimen has an asymp- 
totic cumulative distribution function of the form 
q 1] - exp - , (b > 0), (7) 
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which has come to be known as''The Weibull Function", due tothe em- 


piricalanalyses which W Weibull, of Sweden, carriedouton the strength 
and life of materials by using this function as a basis. Assumption! can 
be justified completely from (7) simply by solving it for N ;. We have 
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*Freudenthal, A. M. & Gumbel, E. J.: On the Statistical Interpretation 
of Fatigue Tests, Proceedings of the RoyalSociety, A, Vol. 216, p. 309, 
(1953). 
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The case of zero endurance limit and zero minimum life. 


In case N, < = 0 and L = 0, the S-N equation (6) reduces to 
,\m 
N..*N, yi— , 
q,S5 qg,U S 
which means thata plotof log N, .< vs. log S willbe linearand ofthe same 
slope for every probability level gq. This is what happens in the case of 
ball bearing fatigue. 








(b) A formula for minimum life in terms of stressand characteristic life. 


Since (7) holds for all values of S, it follows that 
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EXPANSIONS IN SERIES OF 
NON-OR THOGONAL EIGENFUNCTIONS 


Robert Davies 
Research Laboratories Division, General Motors Corporation 


Detroit, Michigan 


ABSTRACT 


Usually eigenfunctions are orthogonal, and the orthogonality is re- 
quired in order to obtain, by the customary schemes, an expansion ina 
series of eigenfunctions. For example, the orthogonality is used to expand 
in a Fourier series or ina series of Bessel, Legendre, or Tschebyscheff 
functions. However, some physical problems lead to eigenfunctions which 
are not orthogonal, and other methods must be used to obtain the eigen- 
function expansion. The paper discusses some of the techniques which 
are available. 


INTRODUCTION 


In many cases it is possibleto reduce a second-order partial-differ- 
ential-equation boundary -value problem toanordinary differential equation 
system by separating the variables. 1-4 Usually a Sturm-Liouville 
system, 5 


= K(x)y! (x), + ‘Ag (x) - L(x) | y(x) = 0, 


a) y(a) + a),2y'(a) + b;,y(b) + b; 2y'(b) = @, 


a. 


) 


iy(a) + aooy'(a) + bo )y(b) + b»oy'(b) 0 


is obtained. In the equation, k(x), g(x), 2(x), the a's, andthe b's are 
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known and Ais an arbitrary constant. The function y(x)=0 is a trivial 
solution forall valuesof A. There will be non-trivial solutions, the eigen- 
functions, for only certain values of A, the eigenvalues. 
If 
aj] 2)]2 by) Pe 
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the system is calledself-adjoint. For twoeigenfunctions, yj(x) and y;(x), 
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Consequently, for self-adjoint systems, 








If A; and <A; are the eigenvalues corresponding tothetwo solutions 
y ; (x) and yj(x), then, integrating by parts, 


b b b 
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y:(a) y ja) | ¥i0) y jb) 
In case A; att the integral itselfis zero. If there is only one eigen- 
function for each eigenvalue, then in case Aj = dj. the functions y ; (x) and 
y j(x) are equal. The integral ~,b 


g(x)y ;(x)dx is positive if g(x)is positive 








in the interval (a,b); then each eigenfunction y;(x) can be normalized by 
dividing by the square rootof theintegral. For normalized eigenfunctions, 
~b 
g(x)y;(x)y j(x)dx = 65; (1) 
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a 
in which ij equals 1 if i = j; otherwise FF is 0 


Equation (1) establishes the orthogonality of the eigenfunctions. Rep- 
resenting the integral -b 
g(x)y; (x)y j(x)dx as aninner product (y;, yj) the 


a 
equation (1) becomes 


(yi. yj) = Oaj- 


If it is desired to have an eigenfunction series 
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xX 7 
f(x) = D_ fiyi (x) 
4=Q9 


— 
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that converges to an arbitrary function f(x), then one can find the coef- 
ficients f; by assuming thatthe series converges and that term-wise inte- 
gration is justified: 


b _b foe) CO b 
f —_ 
[ elxdt(xdy ;(xdax =/ gdy je) fyyGdax = 37 & / gldyiGdyj(xdax - 
Jy “e i=0 i=0 Ya 
fee) 
= >, fi bij; = f;. 
i=0 . ‘ 


Consequently, the coefficients in the expansion (2) should be 


b 
f= | |(x)i(x)y;(x)ax. (3) 
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In inner-product notation, 


f; = (f, y;)- 


ll 








It has been proved that an expansion in eigenfunctions of a Sturm- 
Liouville system of the type athand converges to f(x) ontheinterval a<x<b, 
if f(x) is piecewise differentiable. © 


As an example of an eigenfunction expansion, the Sturm -Liouville 
system 


y(7T) = 0 


yields the eigenfunctions 


Mm 


nV 7 mh ri n 


sin nx sin mx dx 
0 
The corresponding eigenfunction expansion is the Fourier sine series 


OD ~ 7 


f(x) sin nx dx. 


Expansion of the function in a non-orthogonal series. From time to 





time a physical problem leads to a modified Sturm-Liouville system for 
which the orthogonality condition (1) does not hold. Instead there is a con- 
dition of the type 


J] sbelvilelyjeddx ~ aysladyj(a) + Byi(by5(b) = 855 (4) 
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It is assumed that the expression 


b 


rn 
> 


g(x)yf (x)dx - Ay> (a) + By 
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is positive for allof the eigenfunctions and that they can be normalized by 

dividing by its square root. As before, an inner product can be defined, 
b 

j) > E ety, (xy s(x)ax - Ayj(a)yj(a) + Byj(b)yj(b), 
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(y; o¥ 
- 


(yisyj) = 04; - 


From analogy with the previous case, the coefficients inan expansion 


f(x) = Sty, (x) (5) 
i=0 
should be / 
f; = f(x) sy; (x) ] ' (6) 


To derive the formula (6), assume thatthe series convergesand that term- 
wise integration is justified. 


~b 
g (x) f(x)y j(x)dx ~ Af(a)y ;(a) + Bf(b)y ;(b) = 
a 
~b re) fe) ond 
=| g(xdy (x) > fy; (xddx - Ayj(a)> Gyy(a) + By;(b) > “fii (b) = 
J i-0 i=@ i=0 
a 
co b 
= 576! | eGdys(xdy;(x)ax - Ay;(a)y;(a) + By; (b)y;(b) | + 
1=O0 J 
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fe) 
= > f S = f 
ee ee 
i=9 


Consequently the coefficients for the expansion (5) should be 
.b 
, | , : 
f= | g(x)t(xdy;(x)ax - Af(a)yi(a) + BE(d)y;(b) . 


~ 
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which is the equivalent of the equation (6). The convergence of series of 
non-orthogonal eigenfunctions has been established for onlya few cases. 8-10 
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Expansion of the derivative in anorthogonalseries. Oftena physical 





problem leads to non-orthogonal eigenfunctions because the boundary con- 
ditions contain secondderivatives. Differentiating the differential equation 
once will yielda system with orthogonal eigenfunctionsifthe coefficients 
in the equation are constants and if the undifferentiated function does not 
occur in the boundary conditions. Consequently, the derivatives of the 
original non-orthogonal eigenfunctions are orthogonal. Inother situations, 
a linear combination of the original eigenfunctions and their first deriva- 
tives, or a higher derivative, may be orthogonal. 


An example of orthogonal derivatives is furnished by the longitudinal 
vibration of a rod which is clamped at one end, is attached to a mass at 


the other end, and is otherwise free from externalforces. The differential 
equation is ti 


2 e 
Ow _E Ow 


1t< yx 


in which w(x,t) is the displacement of the rod at location x and timet, 
and E and are constants for a particular rod. The boundary conditions 
are 

w(0,t) = 0, Mwy;(2, t) + Ew,(2, t) = 0, 
in which M isthe attached massand iis the length of the rod. If the dif- 
ferential equation is satisfied at x=), then the second boundary condition 


can be written 


“ 


BD xx (2, t) + wy (0, t) = 0. 


Upon separating the variables, the ordinary differential system be- 


comes 
y''(x) + Ay(x) = 0, 
y(0) = 0, 
Ay" (2) + y'(Q) = 0, 


in which Aisa constant. The system does not have orthogonal eigen- 
functions, but on defining 


z(x) x y'(x), 
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z(x) must be a solution of 


Az'(Q) + z(Q) = 0. 


Since a Sturm-Liouville system is obtained, 


the eigenfunctions are 
orthogonal. Representing them by z;(x), 


J 
J z;(x)zj(x)dx =§ 


“0 


and the expansion associated with a function f'(x) is, from the equation (3), 


ij’ 


ll 
f'(x)~ > fiz. (x), fi = | £'(x)z,(x)dx. 


ee... 


If f'(x) is the derivative of a function f(x) whose expansion in a series 
of the eigenfunctions y;(x) is desired, and if term-wise integrationis justi- 
fied, then 


oO @ ¢: 
f(x) | f'(x)dx = >» f’ | i (x)dx => y; (x). 


The series converges to 0 at x=0 and therefore f(0)=0. The desired coef- 
coefficients are 





f' 
= . Integration by parts gives them in terms of the 
1 
functions y; (x): 
' ! y 

f. lr 1 f y or 
f, =F = fh f'(x)z.(x)dx =; — f(x)z.(x)| - f(x)z!(x)dx : 

1Aj VAj J VAj J 

0 ot 








| £(x)y,(x)dx + Af(D)y,() 


*%; | yi) J £(x)y''(x)dx 
% 





The equation (7) yields the same form. However, if the equation (7) 
is used, convergence of the series can be proved easily only for the few 
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cases treatedinthe literature. The methodusedhere ensures convergence 


in case f'(x) is piecewise differentiable and iftheintegration is justified. 


There are special techniques available for justifying the termwise inte- 


gration of eigenfunction expansions 


Wi 
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ON SOME NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS* 


S.D. Conte - Dept. of Mathematics, Wayne University 


ABSTRACT 


This note presents a method of obtaining exact solutions to certain 
types of quasi-linear partial differential equations. Given a partial differ- 
ential equation L(u) = 0, a transformation u = F(@)is introduced. The 
function 8 is chosen as the solution of a certain linear partial differential 
equation, while the function F satisfies an ordinary non-linear equation. 
The proper choise of F and @ frequently leads toa solution of L(u) = 0. 
Several illustrations are given 


The usual linear theory used to describe phenomena in the fields of 
elasticity, hydrodynamics, aerodynamics, and heat conduction frequently 
fails to yield satisfactory results. In such cases recourse must be made 
to non-linear equations, but here the ordinary methods of obtaining so- 


lutions to linear equations, such as separation of variables and transform 


methods, are no longer effective Indeed, exact solutions of non-linear 
equations are very rare. The purposeof this paper is to presenta method 
of obtaining exact solutions to certain classes of quasi-linear partial differ - 


ential equations. 


If we are dealing with a partial differential equation whose dependent 
variable is u(x,t), the methodin its simplestform consists of introducinga 
transformationu = F |@(x,t)| and choosing @(x, t)as the solutionof a certain 
linear equation The function F(8) will usually satisfy an ordinary non- 
linear equation with 8 as the independent variable, the non-linearity in 
the original partial differential equation being transferred to the function 
F. The proper choice of F and @ will sometimes yield a solution of the 
equation inu. The method can best be explained by illustration 


Let us consider the equation 


ea sku. 
t xx 


+ A(u, u,) (1) 





* Read at the regular meeting of the Society, March 12, 1953. 
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where differentiation is indicated by subscripts. It is assumed that k is 


a constant and Aisa non-linear function of uandu,. This may be in- 
terpreted as a generalized equation of one-dimensional heat flow. Some 
important cases of (1) which have been investigatedare those with A = e¥, 


‘ 


e7t/u, uu,. If A uu,, for example, (1) is the equation of heat conduction 


in a moving medium, or again (1) with A lu, may be interpreted as a 
model equation for turbulent flow in hydrodynamics as J. Burgers!h 


’ | 2 
shown We will consider in more detail equation (1) with A = P(u) uf, 
' x 


where P(u) is any function of u, linearor non-linear. We seeka solution 


of (2) in the form u = F |@(x,t)| , the functions F and ® being at this point 


unspecified Ther 
10 7 
uy i vy . u, i 0, ; 
, ) 
8 HIgs 
R. F'6,, + F “e 


and substitution in (2) yields 


) 


F'@, = k(F'9.. + F'6°) + P(F)F'“e° 
‘ x x 


F' |\6, - ké@ 0° |kF" + P(F)F'*| (3) 
Now equation (3) will be formally satisfied if we take 9 and F to satisfy 


the respective equations 


fe) : é... ) (4) 


Of course, (4) is the linear equation of one dimensional heat flow, and its 
general solution is well known“. Furthermore, (5) is an ordinary non- 
oy ar bo ation whose solutioncan be obtained readil Indeed, if we divide 
(5 F' and k, we obtain 


n> Sees # 0 
F' k ' 


an equation which is separable and whose solution 1s 


c19+c9 g(F*‘ = exp|— P/FidF! dF 




















In the simplest case, P(u) equal to a constant, the equation (1) is of con- 
siderable importance and its solution is given by 


5 


u(x, t) F | 0(x,t)|=k In [c,@+4 co | 


= a ca 


as J. Cole? has shown. It may be worth while to consider one case in 


more detail. For this purpose, consider (2) with P(u) | 





u Equation 
(2) may be written 
uu, K uu, + uy : (6) 
and equation (5) becomes 
kK FF" + F'“ = 0 
The solution of the latter equation is 
_k_ 
7 k+l 
F sth (c,8+ co) (7) 
k ] 2 
the constants c,, c. being arbitrary. Thus if 0(x,t)is a solution of the 


heat equation (4), then u = F(6) with F given by (7) provides formally a 
solution of the quasi-linear equation (6). More explicitly, consider the 


boundary value problem consisting of (6) with k = 1 and the conditions 
u(x, 0) = f(x) , 0< x<] 
(8) 
= 0 
u_ (0, t) u_ (1, t) 0 t 
12, 
where f(x) is a suitably restricted function. Then u [20(x, t)} is a so- 
lution of this boundary value problem, the corresponding boundary value 
problem for (x,t) being 
8, : OL. O< x a 9s Ce . 
Q(x, 0) = | f(x) |“ = g(x) , O<x<l , (9) 
6,(0,t) = O6,(1,t)=0 , t >0 


The system (9) may beinterpreted as the boundary value problem for the 
temperature within an infinite slab whose faces x = 0 and x = 1 are insu- 


lated and whose initial temperature is 6(x, 0) = g(x). It is known that (9) 
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*») 
bounded solution~ which 


possesses a uniquée 


may be written 


OO ( ) ) 
TS sexxp(-n*“ 77 t) 
O(x, t) i @ Ch cos n77x (10) 
ns] 
where 
] ~] 
l > 
Co g(x)dx {“(x)dx , 
0 U 
a} l 
Cy ~ (x) cos n77x dx f“(x) cos nix dx 
0 


solution of (6) i 


u(x, t) \2 Co t+ - exp(-n TY cos nix (11) 


and it may be verifieddirectly that the conditions (8) are satisfied. Since 
yntinuous over the domain R(0<x< 1, t >0), then so are 
x» Uy, provided that 6(x,t) does domain R. The 


seen in differentiating 


not vanish over the 


latter restriction may be 


The uniqueness of the 


t) from that of 6(x,t). A more 
detailed discussion « h } 


solution (11) will not be at- 
tempted here 


Another important family 


equations may be written 


(12) 


where f(y at x) isa 


) ) 
nor near function. Whenf =u(v" +¥"), for example, 
x y 
the function 8 is found to satisfy the two-dimensional Ldplacian equation. 
A solution of (12) is then given by u 


-(1n 6) /4 
The method 


yf attack proposed here may be applied to partial differ- 
ential equations of any order, 
; 


but it is obviously limited to equations wit! 


+} 
certain types of non-linear terms 


When it is not directly applicable, a 


yreliminary transformation will often render the equation tractable to this 
4 
method. 
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ELEMENTS OF ATOM] 


C ENERGY 


ANALYSIS AND SYNTHESIS 


Urner Liddel Bendix Av 


This tutorial paper traces the recent dev: pment of physics from the 
discovery of radioactivity to the atomic bom! A st of the known parti- 
cles indicates that over 15 different ones have now beenidentified or sus- 
pected. 

In discussing the problemof aton é rgy t might be fruitful to use 
an historical approach to the problen which is in effect the application 
of an analytical technique to the problen Let us consider the character- 
istics of matter The first property which was discovered was undoubtedly 
that of mass The second characteristic which was discovered was that 
of electrical charge ind the third was the property of magnetism These 
three tundamenta! nara eristics remained tne s e basic properties ol 
matter until the ery recent development r oO! dering the nucleus. 

Matter exhibited itself in three phas« solid, liquid, and gaseous; 
but most of the properties could be ribed by considering the solid 
phase of matter Newton developed s la n the nservation of mass, 
momentum, et n considerationo!l t t iS¢ By t! Sse concepts matter 
was hard and impenetrable 

In examining the ynductivit ( € canal rays (positive ions) 
were discovere’ ind it further appeared that electrical conduction was ac- 

omplished by negatively harged part é . ch ere called electrons. 

Thomson calculated the radius of these part as 2 x iG 13 centimeters 
by assuming that their rest mass was ntirely duet electromagnetic 
harge 

Discovery of radioactivity by Be ierel in 1896 was the first evidence 
of breakdown in the inv lability of matter Hers ot course it was dis- 

overedthat matter spontaneously ad ymposed in such a fashionthat mass 





* Read at the regular meeting ol he So 








was no longer conserved. Rutherford's classical experiments on the scat- 


tering of radioactive particles by thin gold foils showed that there was 
actually a considerable amountof free space separating the atoms of sup- 
posedly solid matter. This space in fact occupied a major portion of the 


volume of the solid. 


Planckin 1900 discoveredthat light, originally considered to be con- 
tinuous, mustbe discontinuous. Thus, we are confronted witha situation 
where two previous concepts of continuum had to be abandoned for a new 
concept of particulate nature. 


In 1905 Einstein published a paper on his specialtheory of relativity. 
Using the Lorentz-Fitzgerald formula for the contraction of electrons in 
fast motionand Newton's second law of motion, hearrivedathis now famous 
equation, E = MC*. This derivation can be accomplished in the following 


manner: 


dE F - dx, since work = force - distance, (1) 
, d(mv) ,. : , ‘ 
also, fk oe from Newton's second law ; (2) 
at 
therefore, 
d(mv) 
ee (3) 
at 
’ dx 
Since — =v, (3) may be rewritten 
at 


dE = v d(mv) 


> 


or dE v<“ dm + mv dv. (4) 


the Lorentz-Fitzgerald equation for explaining the apparent contraction 


of rapidly moving bodies states 








where m = apparent mass, 


m = rest mass, 


c = velocity of light 








w 


Equation (5) may be rearranged as follows 


) i) ) ) e ) ‘ 
= i a yo - = _ 
m“(« ) m~ ¢ j (6) 


Substituting (4) in (7) we have 


dE c*- dm (8) 

( . c (c 
yr E m< (9) 
If light were particulate in nature, it should show the properties of 
momentum, and this was, indeed, proved These experiments led to the 


necessity of the concept of DUALISM. Thisconcept is with us today; and, 


moreand more, we are forced to describe the manner in which we examine 
a physical system to have confidence in the apparent experimental results. 
Thatis to say, the system itselfappears to be affected by the mere process 


of examination 


A logical consequence of this concept of dualism was that matter in 
motion should have associated with it a particular wave length. This wave 


length was calculated by DeBroglie in 1923 to bs 


where h =Planck's constant, 6.6 x 10°°' ergsec. This theoretical concept 
was teautifully proved by the diffraction experimentsof Davisson and 


Germer with electrons. 


This concept of wave motionof matter immediately led to the develop- 
ment by Schroedinger of what isnow knownas the wave mechanics. By this 


procedure, one can calculate the probability that one will find a particle 


ofmatter ata givenplaceata giventime, ||“, from the general equation: 


where 


Y= f(x,y,z) exp(27Ti W t) 


The last of what might be called the historical concepts was developed 











by Heisenbergin 1927 and continues by the title of the ''Uncertainty Prin- 


ciple. By this concept, itis not possible to define precisely the momentum 
and the coordinates of a particle atthe same instant. If one precisely de- 
fines the momentum, the uncertainty in the coordinates becomes equal to 
h. 


The equation is 
SP Ag=zh 


where Ap and Aq are uncertainties in the momenta and coordinates, re- 


spectively. 


Now we are able to see the profound difference between the macro- 
scopic and microscopic properties of matter. The classical mechanics 
deals entirely with macroscopic properties and is governedentirely by the 
concept of continuous fields. The new mathematics dealing with atomics 
is calledquantum mechanics, which describes the microscopic properties 
of matter and deals entirely with discontinuous fields. A corollary deduc- 
tion from quantum mechanics is itsinability to differentiate between obser- 
vationand fact. Wemustuse the equations of probability to define a system. 


As you know, the structure of atoms is today described by a concept 
originally due to Bohr. He conceivedof atoms as miniature astronomical 
systems where the heavy nuclei contain the major portion ofthe mass of 
the system, and these nucleiare surrounded by orbital electrons in a plane- 
tary manner. In astronomical systems, we find only one planet in a par- 
ticular orbit. Inatomic structure, it is conceived that two electrons occupy 
each orbit. The characteristics of the electrons surrounding an atomic 
nucleus are uniquely described by quantum numbers. According to Pauli's 
exclusion principle, no two electrons in any atomic system may have i- 
dentical quantum numbers. Also, the chemical properties of all atoms 
are uniquely described by the quantum numbers of their electrons. 


If we wish to consider a system as composed of particles, we now 
use what is known as the Fermi-Dirac statistics. Here, the y terms in 
the Schroedinger equation are antisymmetrical on interchange of the parti- 
cles. If the y terms are symmetrical, the system is continuous and may 
be composed only of photons instead of particles, and one then uses what 


is known as the Bose-Einstein statistics. 


Letus turn from the discussionofthe electronic structure of atoms and 
consider now the nucleus, whichcontains nearly all the mass and thus all 


the energy of the atoms. Until the end of the 19th century, atoms were 








conceived to be hard spheres The Bohr model placed essentially all the 


massin the nucleus, and it was considered 
tons andelectrons. 
to hold the protons together 
that there could be no electrons in the nu¢ 

from Beta decay that electrons were eje< 
active atoms. The absence of electrons i1 
++ y 


of a new elementary particle called the 


be a relatively stabl:« 


particle composed of 


in 1920to be composedof pro- 
The electrons provided the electrostatic force required 


Later itwas concludedon theoretical grounds 


us, even though it was known 
d from the nucleus in radio- 
the nucleus led to the concept 
itron which was supposed to 


one proton and one electron. 


These particles were discovered in 1932 by Chadwick from the bombard- 


ment of lithium by Alphaparticles. This c: 


the nature of the forces which held the parti 


The true nature of these forces occupies the 


even at the present tims 


The Japanese physicist Yukawa conc: 
perhaps anadditional particle present in the 


discovered. He t] 


ymputed the mass of 
the uncertainty principle andthe DeBrogli: 
is assumed to be of nuclear dimensions 


The mass of the particle may be calculated? 


masses. 
h 
mm 
_ 
6.6x 107° 
2x 107+? Ps 
6.3xmx 3x 10 
4 . » 7 
m 1760 = 10" gramme 
»-00 electronmasses (m 


Such a particle was discovered in the debr 


ation by Street, Stevenson, Anderson, and Neddermeyer in 1938, and here 


the picture remained static until after the 


Wa Ve engtl 


cept left completely unsolved 
cles together at the nucleus 


attention of nuclear physicists 


d the notion that there was 
nucleus which had notyet been 
new particle on the basis of 


1 


The wave length 


TI 
hatis 2 x 10-13 centimeters. 


ee 


pproximately 200 electron 


is resulting from cosmic radi- 


var 


All matter is today consideredto be composedof elementary particles. 


These particles are, by definition, elements of matter which contain no 


structure. Most of them are to be found 


the interaction of supposedly elementary particles with nuclei 


to consider these, it is necessary to define 


various nuclei or result from 


In order 


certain new terms 
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(a) The rest massofa particle is the mass which the particle would 


lave if its velocity were zero. This term is of importance in nuclear 
reactions since they occur at high relative velocities and the momentum 
changes may be greatly different than if considered solely on the basis of 
interaction at low velocity. 


(b) The mean lifeof a particle is the length of time which the particle 
would live on the average if it had zero velocity. Itis obvious that if the 
particle is traveling with velocities near that of light, the life time may 
be greatly extended. The equation for this is identical in form to the e- 


quation for relative mass 


(c) The cross-section of an elementary particle for a given reaction 
is defined as that area which would theoretically result in collision with 
the same probability as disruptive interaction actually occurs. The geo- 
metrical cross-section of atomic nuclei may be calculated on the basis 
of the following formula: 


O= 1.42 10713 A? cm 


where Ais theatomic weight. Experimental cross-sections almost never 
coincide with the geometrical cross-section, varying from .001 to 10° 
times the geometrical cross-section. 


(d) The spin of electrons is evidenced, for example, in the different 
physical properties of ortho and para hydrogen. In the nucleus, the spins 
of the elementary particles are assumed to be additive or subtractive. 


(e) The isotopic spin formalism is perhaps the most difficult to de- 
scribe and understand. It arose from a consideration of the problem of 
the deuteron. This nucleus is composed of one proton and one neutron. 
Consideration ofthe electrical and magnetic forces between these two parti- 
cles leads to forces insufficient to account for the stability of the deuteron 
nucleus. Thus, it was assumed that neutrons and protons were different 
states of the same nucleon. 


It must be emphasized thatno realimportance can be assigned to this 
isotopic spin formalism. Itis useful in considering the theoretical ap- 
proach to particle interaction studies. One danger of placing too much 
emphasis onthis concept may be illustrated as follows. In 1920 electrons 
were assumed to be in the nucleus to account for Beta decay. In 1950 
neutrons are assumed to be in the nucleus, likewise to explain Beta decay. 


One wonders what future concept of elementary particles in the nucleus 


may be derived to account for this same phenomenon. 





| 








a The number and characterist sot elementary particlesare changing 


r at sucha rapic rate that one nesitates totanvulate them w thout emphasizing 

n thatthis ispurelypr naltabulatior Hardlya gr p meeting of physi- 

»f cists occurs without a report of a new par Tables 1 and 2 of ele- 
mentary parti es represented the best tate of experimental and theo- 
retical knowledge in tl! spring of i935 

Le 

e In the early days of study of nuclear phys 5, one wa constrained to 


LV utilize parti« € rom naturai rad act ' T} yoe ol xperimentation 


- has largely ceased With the development of the high energy particle ac 
sla t re } y + +} m< an ” f{ 41 . ++ . . ah! 
ce.ierator , SU eq u¢ O é f ( ( \ OT we are now adie 
to get particies n ty aboratorv with ener é dD & -.5 billion electron 
yn volts. Except for the neutrons discovered by Chadwick, all other unstabl« 
th elementary particles have been d Oovered in cosmic radiation The low 
)- intensity of cOSmic radia nanatne ditt llty Oo arr ng out experiments 
is in altitudes up to 100, 000 feet. se! mitedthe experimentalinformation 
available from these part With +} rato} rer b] 
avallaDbDle irom ese Dat cle e acceierat ve are now able to 
generate nearly a elementary part mn the iaboratory so that it is 


unlikely that we . have t lepend in the future por osmic radiation 

ior our Kn VLE QD ¢ ementary Tra ¢ ‘ a Or However SiInCE€ 

cosmic ray par Le have ne! es p to m n times that which man 
er can generate I il ; eraior ur - é ] i sud A ast tor someyears 
0° to come 


+7 


| Thereare many wa\ f detecting elemental! particles A of these 


nt ways depend upon the nizing properties f charged particle The only 


ns methodavailable for det« ting uncnarged particle is to ol erve the recoil 
e. of a charged particle ejected onco n withthe uncharged particle The 
1our primary t 5 Mm particle detection are 
e- 
Ui rene ‘ rab! =) + 
n > I it n ¢ r ! 
tle f Zat1or & imper°rs 
on Cloud Chan r 
ant | 
The photographic plate wa the original method of detecting charged 
particles. It was usedby Becquerelin ! liscovery of radioactivity. The 
11S te< hnig ue nas markedly advanced 1n tne ast few years and today remains 
p- the best method of observing characteris? f elementary particles. The 
ch scintillation counter was developed by Rutherfordandis usedinmany forms 
ns today We have progressed far from the zir fide crystals used by 
50 Rutherford and today may construct our detecting rystals'' in any form 
v. desired by incorporatior fa suitable phosphor ina transparent plastic 
us Ionization chambers and the more sp« 1lized Geiger counters are gener 
ally used for detection purposes nly The or nal sud chamber was 
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developed in 1923 by Wilson and has progressed today to the large con- 


tinuously recording diffusion cloud chambers. 


Inorder to show some dramatic examples of interaction of elementary 
particles with matter, three cloud chamber photographs are reproduced 
herein. The first, Figure 1, shows the traversal of four lead plates 1/2 


c 


1 


inch thick by a single elementary particle before it is slowed down to such 
an extentthat it caninteract withthe lead. The largenumber ofionization 
particle tracks shown in the four lower sections illustrate the large number 
yf consequences which canoccur from a single primary interaction. This 
is perhaps more dramatically illustrated in the second picture, Figure 2, 
which shows an enormous electron shower generated in the wall of the 
cloud chamber. The last picture, Figure 3, illustrates a true nuclear 
explosion. Here, a single cosmic ray particle enters the photographic 
emulsion, seen in the upper right-hand corner of the picture. It collides 
with a single atomic nucleus, probably silver or bromine, and some 80 
separate particlesare ejectedinthis one explosion. It has been estimated 


that 10 ergs of energy were released in this one interaction. 


Many nuclear physicists agree that the greatest need in the nuclear 
field at the presenttime is a conceptof the nucleus similar tothat of Bohr 
for the atom. There are three models under discussion at the present 
time with differing degrees of support. The liquid drop model assumes 
that the particles are held together by nuclear forces with free motion 
inhibited by viscous flow of the particles. Thus, fission is said to occur 
whenthe surface tension of the liquid is insufficient to hold the elementary 
particles together any longer. 


The powder modelassumes that the particles may act essentially 
independently of each other and are held together essentially by gravi- 
tational type iorces. 


[The shell model has received perhaps the widest acceptance at the 
present time. Here, it is assumed that the elementary particles arein 
energy states comparable to the electron shells in the outer orbits ofthe 
atom. This model wellaccounts for many energy changes found in nuclear 
interactions. Inaddition, it has been discoveredthatmesons, onapproach- 
ing an atomic nucleus, apparently fall into an energy level on capture by 
the nucleus, without falling directly into the nucleus itself. Subsequent to 
this capture, the mesonthen goes into the nucleus with emissionof charac- 
teristic radiation. 


Philosophically, wearein a dilemma. By the uncertainty principle we 


are not allowed to discuss lengths less than 10-13 centimeters; yet the 
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Figure 1 - High Speed Elementary Particles 
Interacting with Lead Plates 
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dimensions involved in interaction of particles in the nucleus must occur 


at dimensions considerably less thanthis. It appears that we shall have to 
abandon the concept of particles completely within these small dimensions 


and say that they have reality only under isolated conditions. 


One reads muchin the newspapers concerning the utilization of nuclear 
energy. The energy involved in these discussions is that of the fission of 
Uranium 235. It can be shown thatif one could obtain the complete fission 
of one gram of U-235, one would obtain an energy equivalent to 2.3 x 104 
kilowatt hours. This is approximately equal to one million watts for one 
day. Itis, of course, extremely difficult to obtain the complete fission 
of a single gram of material. One must, in fact, have many kilograms 
in order to obtain this quantity of energy. It should also be pointed out 
that this energy is essentially thermal energy andis notdirectly convertible 
into electrical energy. On the other hand, it may be interesting to point 
out that one obtains the order of .01 kilowatt hours per gram of coal, so 


that uranium provides nearly one million times as much energy as coal. 


For the last 30 years, physicists have dreamedof the hydrogen conden- 
sation reaction in which four hydrogen atoms coalesce into one helium 
atom. This is believed by the astronomers to be the primary source of 
energy of the stars. If this reaction could be made to occur, one could 
ybtain 1.25 x 10° kilowatt hours per gram of hydrogen. This condensation 
should have the possibility of being operatedas a controlledreaction, free 
from the dangers of explosionand radiation hazards developed in uranium 
reactors. It is perhaps more than a fond hope that we may See this reaction 


as a source of energy within our lifetime. 


In order to finish this cycle, it is perhaps wellto note that ifwe could 
completely disintegrate matter into energy we would obtainfrom one gram 
of any kind of matter 2.5 x 10’ kilowatt hours. It is obvious from these 
figures that any of the reactions which we are able to control utilize only 
a small fraction of the total energy available. Whether we are ever able 
to utilize total energy available from mass remains tobe seen. The mesons 
previously discussed may holda key notonly to the forces holding the nucle- 
us together but to a possible mechanism for the complete decomposition 
of matter into energy. Certainly, if we are ever able to utilize either of 
the last two methods of producing energy, we will no longer be dependent 
on fossil fuels for heat sources. 

In summary then, I hope that it has been shown that matter is rather 
complicated -- the more complicated by our own inability precisely to 
measure the physical quantities invoived. At the moment, the ephemeral 


particles called mesons haveonly significance for the theoretical and ex- 


perimental physicists as items of great interest. However, one cannot 














resist a strong belief that these particles will be effective in the power 
sources of the future. 
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THE INFLUENCE OF AUTOMATIC COMPUTING MACHINES 
ON MATHEMATICAL RESEARCH* 


H.D. Huskey Mathematics Department, Wayne University 


ABSTRACT 


The paper discusses the design and use of high-speed computers, 
and their influence on pure andappliedmathematics. New computer design 
features aimedatmore efficientuse are discussed; anindication of future 
trends is given. Examples in number theory, matrix problems, linear 
programming, and Monte Carlo techniques are mentioned. 


INTRODUCTION 


The influence of automatic computers onmathematical research con- 
tains threeaspects: (1) the opening of a new field in mathematics, namely 
the design anduse of high-speedcomputers, (2) the influence of computers 
on trendsin pure mathematics, and(3)the influence of computers in popu- 
larizing applied mathematics 


Presentand prospective computing machines, have revitalized older 
studies and opened new vistas in many branches of mathematics. A field 
centering around the efficient use of existing computers and the design of 
improved machines has come into being. The development of machinery 
which will precisely carry out complex and intricate logical processes of 
millions uponmillions of steps has ledto the hailing of the eraof the 


' me- 
chanical brains'"' by one school, and to the vehement assertion ofhuman 
‘'corner'' onthe ability tothink by another. In pointof fact, it is exceeding- 
ly difficult to define the word ''think"' in such a manner that an automatic 
computing machine cannot qualify. However, all usersof such equipment 
know that considerable thinking is necessary by the persons concerned 


} 
i 


before these machines can be properly utilized 


The adventof high-speed automatic computers has presented mathe- 


maticians with a challenge to use these machines as efficiently as possible. 





* Read at the regular meeting of the Society, Nov. 14, 1952. 
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The use of computing machines may be considered from two viewpoints: 
programming techniques andcomputer facilities. In any general purpose 
digital computer it is possible so to program that the machine will behave 


like any other specified digital computer. In fact, it has been shown that 
a device with an infinite memory but only a rudimentary arithmetic unit 
can be causedtoimitate any other computer. This means that the signifi- 
cant questionis not whether a given computer can do a given problem but 
merely whether the computer is efficienton that problem. This efficiency 
depends partly upon what special facilities have been builtinto the computer 
and partly upon what programming techniques have been developed by the 
people who use the computer. Thus, this question of efficiency has to do 
both with programming techniques and with actual facilities built into the 
computing device. 


COMPUTER DESIGN 


First consider the computer facilities of the machines now in oper- 
ation. In the initial machines the emphasis was on minimizing tubes and 
simplifying designin order to facilitate construction, checkout, andinstal- 
lation. This ledto computers having a minimum of commands, or atleast 
a minimum variety of commands. One example of this is the National 
Bureau of Standards Western Automatic Computer, SWAC, in which there 
are only eight basic commands. 


Another direction of developmenthas been tolet each digit inthe com- 
mand represent a particular elemental function. A combination of these 
functions thenrepresents anarithmetic operation. Other computers have 
incorporated the arithmetic operations in the memory switching circuits. 
In these cases, a primary motivation was simplicityin design. There are 
computers now inoperation using this principle, as for example, the ACE 
in England. 


With the actual construction of automatic computers and their use in 
computational problems, there have appeared recently machines of more 
complicated design than these first models. This is partly because a 
major bottleneck in preparing problems for solution on high-speed com- 
puters is that ofprogramming. Programming requires trained personnel 
and even then is a lengthy process, partly because of inherent complexity 
and partly because of the difficulty of detecting and removing errors. In 
the future it is hoped that machines will be designed which will permit 
simpler programming. One method would be to provide facilities for ac- 
complishing substantially more erudite commands than are designed in the 
present line of computers. For example, vector operations could be ac- 
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complished with single commands. One machine has been recently de- 


} 


signed to be especially useful on matrix problem 


At this point itmay be ofinterest to note a few particular instructions, 
Originatedto facilitate coding. The repeat commandin the ERA-1101 com- 
puter is designed to expedite obeying a group of commands a specified 


number of time 


The" B" register usedinthe Ferranti computer designedat Manchester 


University, England, allowsa givensetof commands to beappliedto vari- 
ous sets of number without modification of the ommands. Another aid 
to coding is thx r itive address'' system in the National Physical Labo- 


ratory's ACE (pilot model). Here parts of the memory can be precessed 
relative to other part: o that a given set of commands can be applied 
to sequences Oi number wit only the location of the commands being 


changed. 


The recordcommand of the UNIVAC and the subroutinetransfer com- 
mands of the CIRCLE computer facilitate the entry into subroutines from 
the main problem routin« The problem is how does the computer return 
to the proper place in the main routineafter doinga subprogram operation 
such as a sine or c ne The record command, or subroutine transfer 
Specilies this inor vord, whereasinother computers two words may be 


required 


A more recently designed computer the DYSEAC at the National 
Bureau of Standard ncorporates someinteresting commands It has a 


summation command which sums a sequence of numbers in the memory. 


To facilitate floatir point computations a "justify or normalizing com- 
mand shifts the number into standard form. (In a floating point operation 
the digits ofone wordposition in the computerare dividedinto two groups. 


One group represents the mantissa or significant figures of the number, 
and the other the exponent of the bas« For example, ina decimal com- 
puter +123456789 might representthe number +0.1234507 x 1039, the last 
two digits representing the exponent plus 50.) There are two command 
counters either of which may be used by any command. This facilitates 


subroutine entry and return to the main routine A flexible system of stor- 
ing and resetting these command counter contentsis suppliedbya so-called 
'file'' command. A breakpoint file command provides these same facilities 


1 


but with optional use depending upon the manual setting of a switch 


ine 


PROGRAMMING 


Let us next turnour attention to programmingtechniques. Here isa 


subjectthat lies in-between the engineering and design of computers, and 
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the use of them by mathematicians. It seems preferable to classify pro- 


gramming as a mathematical subject since it is a necessary adjunct of 
operating a computing system. The operation of a computing system is 
primarily a mathematical activity, whether it be involved in engineering 
design orin numbertheory. The history anddevelopmentof programming 
techniques is comparableto thatgiven above concerning extra facilities in 


computing machines 


A high-speedautomatic computer is sucha complex device that methods 
of using a given computer continue to improve and develop even without 
any change in the logical or electronic aspect of that computer. Yor ex- 
ample, just recently there have been devised so-called post-mortem check- 
ing routines which process other routines in order to find out what is wrong 
withthem. Work isalso being done on floating or relative address systems 
that simplify the bookkeeping processes involved in laying out a program 
of computation. In these floating address systems the actual instructions 
involvedin carrying out the computationare classifiedinto various groups, 
and thus, in laying out the program it is only essential to make sure that 
the commands ineach group arg laid out in the proper order. Commands 
can be added to each group without upsetting the general arrangements of 
the other groups Programming techniques of this sort have simplified 
greatly the preparation of problems for computation on high-speed ma- 


chinery. 


In this direction mathematicians may look forward to computing ma- 
chines for which a problem can be essentially programmed by typing the 
formula that is to be evaluated. There would be no worries about where 
informationis to be stored, or what logical patterns are to be involved in 
the computation. This sounds quite wonderful, of course, but on the other 
hand, it is doubtfulif any ofthe purely mathematical preparation that must 
precede the running of a problem can be by-passed. For example, some- 
one must decide what mathematical method is to be used, and how many 
Significant digits are to be carried. Once these questions are decided, 
it is a simple matter to write down equations representing the processes 
involved. The equations can then be typed into the computer. Next the 
computer mustinterpret these equations and produce a list of commands, 
which when processed by the computer will give the answers to the desired 
problem. 


Now consider the effect of automatic computer programming on de- 
velopments in logic. One of the questions, of course, is how does one 
write down the logical equations that represent a computer design? This 
operation may involve multiple valued logical functions, and more par- 
ticularly, the equations of a computer involve timeas anindependent vari- 


able. The techniques for writing the equations that representthe response 
; 
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ofa logical circuit are quite involved, and itis rather difficult to write 
such equations and make sure that they are correct. In general, most of 
the designof digital computing machines has up to this time at least been 
done by a trial and error, or brute force method. But perhaps techniques 
will be developed in applying logics to this field so as to derivea straight- 
forward procedure tor arriving at the optimum design of an electronic 


Omputer or similar de Ce 


APPLICATIONS TO PURE MATHEMATICS 


Consider v th eldof pure mathematics. One of the first types of 


problems pers King out computers like to tryis factoring problems 


This is because here one can specify quite large scale problems requiring 
] long rinnir to Ty ; KNOW)! ths. ne er y | , EF > ws le 
é tong runnin ! ] ] y »W t i vs I naavance Or examp Ce, 
assume the mputir machine is instructed to divide 21% by 2 As itis 
known that thea : the computer can be instructedto start sub- 
tracting 2's fron ° and totally how many timesit does this. It will run 


on for quite a period of tims At the end of the process, by looking 


atthe remaind« I noutla | ( one can check thatthe computer proba- 
bly ran correct ! i i e number I peratior 


Worthwhile ap} tior of computing devices have been made in the 
fieldof factorir rpe rime Of particularinterest is the investigation 


of Mersenne prin that has beer irried on as an incidental project at 


the Institute { Numerical Analysis in Los Angeles (One interesting 


point about this roblem is that it represents one case where a non-pro- 
fessional coder prepared a problem, and the code turned out to be correct 
in every detail.) Ti purpose is to study numbers of the form 2P-1, and 


to find out for what valu: f pthese nun rsare prime (Mersenne primes). 
In general, the pr ram is set up to handle numbers for p's as large as 
2293 The limitationon the size of the number involved is the amount of 
room in the high-speed memory to store the number Also for numbers 
much larger than this omputation time becomes a limiting factor. But 
to summarizein brief the results that were obtained, the table of Mersenne 
primes was extended by the addition of five new primes, including the 


largest known prime numbers at the present time 


Theinterest ir lersenne primes arises from their relation to perfect 
numbers, a perfect number being one which is equal to the sum of all its 
exact divisors; e.g 8 1+2+4+7+14 The only known general expression 
producing perfect numbersis due to one of Euclid's theorems which states 
that 2P-1(2P-1) is perfect if 2-1 is prim«s Mersenne alleged in 1644that 
2P_) 


is prime for certain numbers, p, namely, 2, 3, 5, 7, 13, 17, 19, 


ea an . P 4 Ze ‘ > £ , j fs ‘ ? : 
oi. yy - 1Z7 an Ao ~~) losequently i ana 25/ were itounda not to resuit in 
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primes, and also it was discovered that Mersenne hadomitted three num- 


bers: 61, 89, and 107. In general, only in the lastten years have some of 
the final numbers been laboriously checked. For example, Dr. D. H. 
Lehmer, spending some 700 hours ona desk machine, was able to show 
that 2*9/.] was a composite number. The SWAC was able to reach the 
same conclusionin 48 seconds. TheSWAC on its firsttry on this problem, 
picked up two new Mersenne primes, namely, 2521.) and 2697-1. Subse- 
guently, three more numbers were added, the largest being 22281). 


This investigation of large Mersenne primes has other qualities of 
interestto those preparing problems for high-speed automatic computing 
machines. First, note that the largest prime number investigated in this 
system has 2281 digits inits binary representation, and thus, requires 
considerable space to write down. It may be asked, ''How can a computing 
machine with a limited number of registers for storing information be 
used to compute with respect to a number of this size, which requires 
most of the memory just to write down the number?'' The answer illus- 
trates the power of, and the need for, mathematical methods along with 
computing machine methods. In this particular case a theorerm of Lucas 
is used to prove the primality of the numbers involved. Lucas investi- 
gated Mersenne primes, and it is quite likely that he became tired of di- 
viding in order to find outif a particular number was, or was not, prime. 
He eliminated a great deal of the tediousness when he evolved the following 
theorem: To investigate the primality of a number such as 29-1, which 
is equal to 31, construct a Lucas series in which the first term is 4, and 
any other term is the preceding term squaredminus 2, etc. For example, 
if the first term is 4, the second term is 16-2, or 14, the third term is 
190-2, or 194. Actually, one no longer has to deal with 194, since one 
is only interested inthis number relative to 31. Hence, it canbe reduced, 
modulo 31, which means that instead of dealing with 194 one only has to 
deal with the remainder when divided by 3l, namely 8. The 4th Lucas 
numberis 8 squared -2, or 62, etc Now, if 2P-1 is prime the p-1 term 
of the Lucas series is divisible by 2P-1. In this particular example, the 
fourth term is equal to 62 which is divisible by 31; therefore, 31 isaprime 
number. This theorem was usedin checking out the primality ofthe larger 
primes, and the important thing is that the machine had to store only the 
reduced number modulo 2P-1. This example illustrates the power of ad- 
vanced theoretical developments used in conjunction with improved com- 
puting methods and equipment. 


High-speed computing techniques have been applied to other problems, 
for example, the Riemann Zeta hypothesis. No one knows whether this 
hypothesis is true or not - no one has been able to prove it - no one has 
been able to find a counter example. Many numbers have been investi- 


gated, so if counter examples exist, they are known to involve rather large 








_— 








ow 





numbers Certainly one of t} 


chinesis their use to di: ver 
ever, to the author knowled 
chine actually to prove a mati 
has been propose that pos sibl 
simple setof axioms Ould be 
ing all the 1; - t f the 


COMPUTERS IN APPLIED MA 

The fi i mat} 
puting machines, n be divide 
and engineeri! 

In the tal ele + 
tabulation of c¢ is data has ps 
ing machines as m iS al t! 
the OrTriginai ee rnry r ige!l 
computer tor t rt liar 1 
study contract e Ecker 
mately result nm the pr luct 
Bureau This m be 
has bee! iseqa Le I 
fact, the 19 ' 


cessead TY 


Computer to 


Tabula t 


industry ar 


omputir t i —T. 
Latior or rt ’ 
,Ories T} 
yutTput equipm«e j mr) 2a " 
to the ymmpu r <7 
University p soil a 
ment Their pr 
tables and the ENIAC is r 
problem. In fact mpute 
the shell sps« t Witt 
ympute r tne t nmor 
point I é r b } . 


imMpi¢ to such theorems. 


1es Of automatic computing ma- 


How- 
has yet used a computing ma- 
this sort, although it 
lepending upona reasonably 


e computer simply by exhaust- 


onsidered in relation to com- 


statistical 


e Bureau of the Census in the 


redthe deveiopment of comput- 


n that this Bureau was one of 
ras DI ite industry to builda 
re ilt of their interest was a 

n te Corporation which ulti- 
VAC to be ised by the Census 
n for several years now and 
asonabdly large scale In 

no populated states was pro 
ate i months for the Census 

Y I eCOpl_e 

problems that arise in 

est ct o the amount of actual! 
é problems involve tabu- 

( tems in various cCate- 
irpose need high-speed input- 


will supply information 


r, the ENIAC. was builtat the 


-t ofthe Army Ordnance De part- 


od method ot producing tiring 
solution of this 
a shellin less time than 
speed electronic 
trajectory problems from the 


or o represent the motion of 


45 





One of the problems that has been of interest to the aircraft industry 
is the vibration ofairplanes. In general, the study of such vibrations 
involves studying flutter matrices to find their eigenvalues. A great deal 
of workhas beendone at theInstitute for Numerical Analysis to find better 


methods of solving systemsof linear equations, of finding eigenvalues for 
matrices, or for inverting matrices. These methods range from study of 
simple elimination methods to rather elegent methods involving steepest 
descent techniques. Steepest descent methods involve approximating a 
vector solution to the system and modifying this vectorin a direction which 
minimizes a residual, or the sum of the squares of the residuals, when 
this vector is substituted in the system of equations. 


A new branch in appliedmathematics that has developed hand in hand 
with the developmentof high-speed automatic computation is thatof linear 
programming. Theideasandconcepts of linear programming are notnew, 
but in general the problems that are involved in this field have so many 
equations, or inequalities, and so many unknowns, that it was impractical 
to work out examples before the advent of high-speed electronic compu- 
tation. 


With the advent of high speed computing, so-called Monte Carlo 
methods became almost a fad among mathematicians. There is a corre- 
spondence between certain sequences of random events and certain types 
of differentialequations. Thus, the solutionto certain problems involving 
sequences of random events can be found by solving the corresponding 
differential equation. Conversely, certain differential equations can be 
solved by computing sequences of random numbers and processing them 
in a prescribed way. 


This can be illustrated by a problem called the Gambler's Ruin. 
I'wo people G and H play a game involving turns. The probability that G 
wins on any turnis aconstant p, and the probability that H wins isq. The 
stake at each turn is h. What is the probability that G ruins H if their 
combined fortune is b? 


Let x(t) be the fortune of G at any time and let g,,,(x) be the proba- 
bility that G ruins Hin not more than m turns. 


Consider twocases: (1) G wins next turnand ruins H in notmore than 
m-1 turns, (2) G loses next turnand ruins Hin not more than m-1 turns. 
Therefore 


g(x) Pg, -) (xth) + G8 m-) *-h), O0<x<b, g,,(0) = 0, g (>) = i. 


By letting moo and h+0 a differential equation of the form 
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where lim }{ 
obtained 


The fractionof the times that G ruins H in these games may give a solution 


to the differential equation 

Number vith sufficiently random properties may be generated ina 
high-speed computer, and in effect, cause the computer to play a number 
of games such as ds ribed above 


Similar''random walk" techniques cansolve La Place's equationat a 


given point in a region with prescribed boundary values. Most other 
methods of solvi: ch equations give answers over the whole region, 
and may involve n hmore computing than that needed if the value of the 
solution at one point ill that is desired 

In closing, I find the following points are worth emphasizing. First 


there are many extremely interesting problems which by presently known 


methods would require years of computing even on the fastest electronic 


computing equipmer ivailable For these problems new methods are 
needed. The 1 point that the process of automatic computation 
is almost ar irt it the present time The time used toprepare and 
check problen mpetes with the time required to run them. Improve- 
ments both in computers and in techniques of using these computers are 


This leads t i final point thatalthough numerous commercial compa- 


nies aredeliveringa matic computing equipment to alltypes of custom- 
ers all over the yuntry there is very little provision for training people 
to use these mac! efficiently However, Wayne University is pio- 


neering in the development of a comprehensive training program, and it 


s the author's hope that in the next few years there will be a substantial 


increase in tnis VI { training activity 
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LINEAR ANALYSIS OF FEEDBACK-CONTROL SYSTEMS* 
E Calvin John On 


Bendix Aviation Corporation, Research Laboratories 


ABSTRACT 
his tutorial paper presents the mathematical background underlying 


the analy SiS Ol teeaback Systems Ope rating i1n the linear region. Classical 
and transform methods of solving the equations are taken upand some dis- 
cussion devoted t stead tate frequency response and its relation to 


system stability 


1. NONLINEARITY 


As the subjectindicats this discussion is confined strictly to linear 
systems. Actually, all physical systems, if pushed hard enough or ex- 
amined closely enough, will exhibit some nonlinearity. Consequently, 
the assumption of linearity is always an approximation. For many pur- 
poses, however, the results so obtainedare quite useful. In anycase, the 
relative ease with whic! near systems can be handled mathematically 
provides considera! ncentive for the use of linear descriptions when 
they areat all appror ate That only the mathematical modeland not the 
actual physical system is linear, however, should not be forgotten. 


2. SYSTEM ELEMENTS 


Most linear elements found commonly in feedback-control systems 
can be characterized by ordinary linear differential equations having real 
onstant coefficients Thus f the box of Figure 1 represents an element 
2i some sort é tric motor 





Figure | System Element 





*Read at meeting of Study Group on Feedback Control Systems Nov. 30, 
1951 
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hydrauli< 


motor, mechanical or 
Output x, 


electrical network, 
as 


a function of time t, often can be related to its input xj(t) by 
a differential equation of the form 


dx 
; i ‘ie 
ie T AoX = bm ¥ * by 
atm 


dt 


n which the coefficients aj and b; are real and constant. 


3. EXAMPLE OF MOTOR 


As an example, consider the d.c. electric motor represented sche- 
matically in Figure 2. 


The motor is driven by applying a 


a e 





Figure 2 - Armature-Controlled d.c. Motor 


. 
voltage 


1eld 


e 


; to the armature terminals 
be 


The field current I¢ is assumed to 
constant so thatthe torque developed by the motor is roughly pro- 
portional to the armature current: 


ki, 


Neglecting 


— 
™~m 
~— 


inductance 


»f the armature circuit, ig is given by 


(3) 
in which Rg is the 


motor The 


armature resistance, and e,, is the back emf ofthe 
back emf in turn is proportional to motor speed 








or whatever -- the 





.e 
'y 


the 





Cm = KX, (4) 


The torque developed by the motor is absorbed by the load, which is as- 
sumed to be purely inertial, J being the inertia. These relations lead to 
the differential equatior 


adWo k‘ k 
oe nto Yo - ej (5) 
dt Re R. 


Equations of this type frequently are normalized to make them more 
generally applicable by defining the new constants T (called the ''time 


constant’'') and K 


3 nae (6) 
K =~ (7) 


The equation then may be writter 





T 2 + w, =Ke (8) 


4. RESPONSE OF MOTOR TO STEP FUNCTION 


If the form of the input voltage e; is specified as a function of time 
the resulting speed wu), canbe found by solving the differential equation (8). 


Suppose, for example that e; has the form shown in Figure 3, i.e., a 
"step function'' of magnitude V. For t >0, (8) becomes 
. 
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ro + KV (9) 


The solution of linear differential equations having constant coefficients, 
such as (9), is the sum of terms of the form eP* which satisfy the homo- 
eneous equation(right member of the equation set equal to zero) and a 
particular solution of the inhomogeneous equation. 


icular solution is the steady-state term KV. The homo- 


. a, 2 0 (10) 


Assume 


in which C is any constant Substitution of (11) into (10) gives 


‘(Tp + 1) C eP* = 0 (13) 


Since eF* is never zero, (1l)1is a solution for nonzero C only if 


ee (15) 


Ci + KV (16) 


if , 1s initially zero, C must be -KV, so that in this case 











O- 


3) 








The form of this solution is sketched in Figure 3. 


5. TRANSFER FUNCTIONS 





d 
In terms of differential operator D = at’ (8) can be written 

(TD + 1)W, = Ke; (18) 

or 
a. (19) 

ej TD + | 
The function 

H (D) (20) 

TD+ ! 


is called the ''transfer function'' of the element. The transfer functionof 


the more general element, described by (1), similarly is 
B (D) 
H (D) _— 2] 
ee) =" (D) (21) 


in which 


A(D) aD + ay D+ a, (22) 
B(D) = b,,D™ + --- + bj D+ b. (23) 


6. SOLUTION OF THE GENERAL EQUATION 
If the input x.(t) to the elementof Figure 1 is given, the right side of 
(1) is a known function of time. Here again the output xo(t) is the sum of 
a particular integral together with terms of the form ePt which satisfy 
A(D) x, = 0 (24) 


With eP* substituted ir (24) for x,, the result is 


A(p) ePt = 0 (25) 
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from which 


A(p) = 0 (26) 


IfAisann th degree polynomial thereare n values of p, p,, which satis- 


y (26), called the characteristic equation. If the Pp, are all distinct, the 


general solution of the homogeneous equation (24) is 


t 
i 


- ' C. ePkt 


n arbitrary constants. 
ROOTS OF THE CHARACTERISTIC EQUATION 


yeneral are complex, say 


(28) 


sin A(D) are real, the py, always occur in conju- 


resulting solution can be real -- of the form 


(29) 


ermine to a large extent the nature of the behavior 


The values of p, det 
of the element In (29), for instance, if any QO, > 0, anoscillatory mode 
is possible in which the amplitude of oscillation increases exponentially 
th tims On the other hand, if 0; <0, theamplitude decreases expo- 
nentially with time For dg k 0, the amplitude remains constant. When 
k 0, the mode is not oscillatory but merely a simple exponential, in- 
reasil yr decreasing depending on whether Q;,>0 ora;,<0O. 
In Figure 4, some possible values of pK are spotted in the complex 
p-plane The behavior of the modes corresponding to the various P, 1s 
sketched nearby. Thus, a certaintype of behavioris associated with each 











Figure 4 - Behavior Associated 
g 


8. STABILITY 


A rigorous, completely sati 


complicated and depend 


SOmewnhat 


With Various Portions of the p-plane 


sfactory definitionof 


stability ''is rather 


+ 


upon the situation. The general idea, 


however, is thatina ab -lement or system the effect of disturbances 


eventually di out Evidently an 


element having excitable modes repre- 


sented by root nthe right half p-plane or onthe imaginary axis would be 
unstable in this s: é 
Much emphasi placed on the question of stability in the analysis 


and designof fee ick-control systems since anunstable system ordinarily 
is completely unsatisfactory. In fact simply to know whether a system is 
stable or not is not usually enough -- it is also important to know how 
stable a devic« More is saidaboutabsolute and relative stability later. 


9. OPERATIONAL OR TRANSFORM METHODS 


The response of the motor 


} } 


to a step function is obtained in Section 4 


by the so-called classical solution of the differential equation. Operational 


methods based on t} Fourier 


tensively in ser\ analysis 


or Laplace transform also are used ex- 
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The Laplace transform of a time function x(t) often is written and 


defined as follows: 
co 
ms x(t) | X(s) x(t)e7S* dt (30) 


O 


Thus, a functionof time x(t) is transformed to a function X(s) of a complex 
variable by the integral of (30). The time function can be recovered by 
means of the corresponding inverse transform: 


-~-A +10 


mw -t | X(s) x(t) X(s) eSt ds (31) 
It follows from these definitions that the transform of (1), assuming 
the time functions and their first few derivatives are zero initially, is 
n ‘ m , 
\a, s** + + ajS + aq) X(s) (b,, Ss Hcoee + b)s + bo) X;(s) 


(32) 


A(s) Xo(s) = B(s) Xj(s) (33) 


In the above, X;(s) and X,(s) are the Laplace transform of the inputand 
putput respectively If xi(t) is known, its transform X;(s) can be found 
either from the integral definition (30) or, as usually is the case, simply 


from previously preparedtables which list various pairs of corresponding 


t and s tunctions Then the transform of the output is 
B(s) : , 
Xo(s) = —S5) x;(s) = H(s) X;(s) (34) 
A(s) 
from which the output as a function of time can be found by taking thein- 


x(t) = 27) |x5(s)| = £7! |H(s) x, (s) (35) 

If the initial values had not all been zero, another term would have 
uppeared on the right side of (32) and (33), which would have been taken 
into account in evaluating the inverse transform. In this lies one of the 
chief differences between the solution of equations by Laplace transform 


methods and by classical methods. The classical approach yields a general 








ga —— 











solution with n arbitrary constants the Cy of (27) -- which must be 


evaluated from tl tial conditions. In the Laplace transform solution, 


on the other hand nitial conditions are taken into account in the trans- 


formation of the differential equation. The resulting time function eventu- 
ally obtained through the inverse transformation, then, is already adjusted 
for the proper initia! conditions 


Both methods of soluti nrequire factoring the characteristic equation 


to find the n p-roots and bot! of courses require a knowledge of the 


initial conditions rhe assical method, as already mentioned, requires 
the evaluatior nstants, whereas the transform method often requires 
a partial-fractior xpat r r the equivalent The transform method 
seems to be some it more systemat ind may provide a neater way of 
obtaining the part ir integral for me input forcing functions. Other 
reasons for it Kter eu ire that some useful results of a general 
nature ar: 311 taines his iy and that transform methods aid in 
unifying transient and uusOidal steady-state concepts. For the actual 
computation : solut to a specif equation, however, the classical 
method ofte1 n e direct and just as tractable as transform methods, 


particularly fo rder yuation 


The trar f in elem«e frequently is written as a function 


of s andcons red e a Laplacetransform Thus, the transfer function 
of the electric n r ire may be writte 


His) (36) 
_ 
[Compare 


10. THE SINUSOIDAL STEADY STAT 


Section 4 r : the esult of applying a step function command 
to the electri rr t | 5 ty f input 1s usea Ire guently as a test sig- 
nal to provide r tion of the performance of an element or system. 
Another standa t nal is the inusoid Suppose, for example, that 
the motor input e;| sinusoidal! say 


e.(t) = E t=Re E, e* (37) 


Because ofthe natu I the diilerentiai equatior (8) governing the motor, 


the steady state respons also is sin [It may be written as 


(t) =Refl. e™ (38) 
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in which £1). may be complex. Substitution of this assumed form of the 


solution into the differential equation (8) gives 

Re(Ti WM, Quis « el ase) - K Re E; el Wt (39) 
Since this must hold for any t, it reduces to 

iw T + 1) 22, KE (40) 


or 


Thus the ratio ofthe complex amplitudes of output and input is simply the 
transfer function, now writtenasafunctionof iu). In this form the transfer 


function sometimes is calledthe sinusoidal or frequency response function. 
ll. FEEDBACK-CONTROL SYSTEMS 


Control systems usually are composed of several elements, each 
having its own transfer function. A typical servo might be represented 
schematicallyas shownin Figure 5. In additionto the differencing element 
the system con- 


which forms the error 6, fromthe input 6; and output 0,, 
€ 1 Oo 


tains three elements with the transfer function Hj(s), H.(s), and H(s). 





Phe r-alltransfer tunction of the forward part 
) Q, 0, 
A 
> > H , (s) - H3(s) ° > 
6. 
Oo 
A 
oo) - € = 
€ 9; 9, 
Figure 5 - Schematic Representation of a Servo 
fi the loop, relating 9, to 6, is the product of the individual transfer 
functions 








H 





—————— (s) is S : is 
; i H,(s) H,(s) H3(s) = H,(s) (41) 


The transform of the error equation is 


@ .(s) = 8;(s) - ®0(s) (42) 


The transfer function which describes the behavior of the complete system 
is obtained from (41) and (42) by eliminating 8 (s) 
e\ 





—— (43) 


In Sections 7 and 8 we saw that the stability of an element depends 
upon the position in the complex plane of the zeros of the denominator of 
the transfer function For the system of Figure 5, these zeros are the 
roots of 


1+ H.(s Q (44) 


H,(s) (45) 


Therefore, inorder that the complete system be stable, the roots of (44) or 


(45) must lie in the left half of the complex plane. 


f the characteristic equation can be factored andthe actual roots ob- 


tained, the stability question is easily settled Usually this would be a 


long andtedious job if the characteristic equation is a polynomial of very 
high degree. And sincethe relation between the polynomial zeros and the 


individual system elements is somewhat indirect, a complete redetermi- 
nation of the zeros may be necessary even for slight modifications of the 


system 


A simpler method might be application of a test such as Routh or 
Hurwitz, whichdetermines the signs of the real parts of the zeros without 


— 


‘ 


the necessity for factoringa highorder polynomial. These tests, however, 
give directly no measure of the relative stability of the system. And here 
again, slight modificationof the system may require completely new calcu- 


lations. 
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13. THE NYQUIST CRITERION 


By far the most useful tool for the design of feedback systems at the 


present time is the Nyquist stability criterion with its variations and ex- 


tensions Use of this criterion involves plotting the open-loop transfer 
function H,(s) in the complex planeas s traces out a closed contour which 
includes the righthalf s-plane. This contourin the s-plane usually has the 


form sketched in Figure 6a which, for large enough R and small enough r, 
enclose all of the righthalf-plane that is of interest. As s traverses the 
contour of Figure 6a, H,(s) takes on various values which when plotted 
produce a closed contour or locus in the Ho(s) plane. A typical locus is 
sketched in Figure 6b. Corresponding values along the s and H, loci are 
labeled Thus, the positive imaginary axis in the s-plane, from large 

to small u leads to the portion of the H, locus fromatob. The small 
semicircular detour around the origin in the s-plane is taken to avoid the 


point s 0 at which the H,(s) for most servo systems becomes infinite. 


The Nyquist criterion states that, if the individual elements whose 
transfer functions make up H,(s) are themselves stable, the complete 
system is stable provided the H,(s) locus does not enclose the -1 point. 
If, for example, the scaleof the sketch of Figure 6b is such that the point 


x corresponds to -l, the Nyquist criterion indicates that the system is 


unstable However, if y corresponds to -l, the system is stable. 
1 4 iln,H, 
i _, , 
—~ 
¥ *. 
v R P \ 
r \ 
} y | a 
= > + a ~ KTS mn ~ ~ ~ > 
. - —“b ReH 
v / 
< 
ad — dtah>-" 
(a) s plane (b) H,(s) plane 


Figure 6 
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As a result of the fact that the coefficients in differential equations 


describing the system are real, the Hp locus forWnegative is always the 
mirror image about the real axis of the locus for positive uw Consequently, 
to apply the Nyquist criterion to the usual servo, it is actually necessary 


to plot only the transfer function H,(iw) aswvaries from very large to 
very small positive values. For stability, then the -1 pointmust lie tothe 
left of the locus. 


14. RELATIVE STABILITY 


One of the greatadvantages of this approachis thatan idea of the rela- 
tive stability of a system can be obtained simply by noticing how closely 
the H(i ) ) locus approaches the -1 point. For reasonable system stability, 


itturns out that the locus should not enter a circle in roughly the position 


shownin Figure 7. A locus plot such as this, then, gives a clue as to how 
the elements of the loop should be designed so that whenthe loop is closed 
the system will perform satisfactorily A little experience shows, for 
example, what modifications or additions might be made to a system pos- 


sessing locus a, Figure 7, so as to produce the more acceptable locus b. 


w 


” Relative 
stability circle 
& 
oll _ 
a b 
¥ 
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15. THE LOGARITHMIC TRANSFORMATION 


A variation in this procedure that is of considerable practical value 
consists of plotting the logarithm of H,(iw) rather than the functionitself. 
The log of H(i J) 1S 


In H,(iw) = ln |H,(iw|+ i Ang H, (io) 


It is customary to plot the log of the magnitude along the vertical axis in 
units of 20 log)g |H,, (3 )| (called decibels whether associated directly with 
power or not) and the angle along the horizontalaxis, as inFigure 8. The 
point 0 decibels, -180° corresponds to the -1 point of Figure 7, and the 
relative stability circle there transforms to the oval sketchedin Figure 8. 


A typical locus also is sketched in Figure 8. 


This transformation simplifies design procedures primarily because 


the magnitudes of the individual transfer functions can be 


Typical locus 


decibels 
oo 
Critical point 
0 eo t 
7 
Transformation of relative 
stability circle 
y 
- 180° -90° 0° 
Figure 8 


multipliedto obtainthe over-alltransfer function H, simply by adding their 
logarithms. Also, there exist simple graphical constructions for loga- 


rithms of the magnitudes of the individual transfer functions. 
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DIREC TIONAL CONTROL OF AUTOMOBILES* 


Motors Corporation 





[The paper describes theoretical and experimental approaches to the 
problem ofdirectional stability of automobiles. For a specific passenger 
ehicle, tl ywest mode of oscillation has calculated frequency in the 
yrder of ps Damping is about .5 of critical under conditions of fixed 
YT ind about -5 under free yntrol These values are confirmed 
ipproximately by experiment. Frequency and damping of the second mode 
ire only roughly predicted by calculation. 

Whenever men have tried to construct a new type of vehicle such as 
the irplane yr the automobile, it has been quite clear to them from the 
be ning that some kind of directional controlis necessary, butonly gradu- 
i y ive tney reaiizedtnat yntrollability alone is not enough; the vehicle 
mus als< | ssess sometni y generally called'' directional stability. By 
tl ern 5 usual meant the ability to follow a desired path without di- 
rectional oscillations Efforts to produce such a property by cut and try 
metnods a ne nave yiten Deen unsatisitactory and some understanding Oo! 
the indamenta lynamics yfter ymbined with considerable analytical 

Tr na pee! I essary 

The ntroductior f softer tires and softer suspensions of the inde- 
per nt typ n Amer in passenger ars led to such radical changes in 
t] r tional properties that research into the fundamentals of steering 

is begu! Valuable information was developed, but evennow the problem 
{ no means fully understood 

The directional yntrol forcesof a vehicle come either from internal 
etfects Ssucn as inertia or gyroscopic torces or from the interaction of 
th el e and the medium in which it travels, such as hydraulic forces 
ma vesse 1erodynan forces on anaircraft, and ground contact forces 
ma road vel ‘ These round ntact forces of pneumatic tires were 

é ited y Evans?‘ and by Bull‘ Their effect on the car's action in 
ir equ riun ynditions was discussed by Olley?: 4and was measured 
Read at the regular meeting of the society, Feb. 12, 1953 











be } ] } } ! as } 5 
by the so-called ''skidpad test'' at the General Motors Proving Ground. ” 


This test furnishes a great deal of information, but since it is a static or 


equilibrium test, it cannot measure the dynamic properties of the car 
which appear in transition maneuvers A dynamic transition test, the 
checkerboard test,'' was also developed, using the rather cumbersome 


instrumentation techniques available at the time 


Recorded data from one such test are shownin Figure 1, whichis a 
reproduction of Figure 19 in reference 5 This lustration shows at the 
bottom the path of the car over the checkerboard, proceeding from left to 


rightand nextabove this, the steering angles measuredat the front wheels 


and at the steering whee! It can be seen that the maneuver consisted in 
subjecting the car toa step function input at the steering wheel, first to 
the rightand then to the left. The upper traces record slip and roll angle 
and finally pathangles and lateralacceleration allas functions of distance 


traveled Note thatit requiredapproximately 60 ftoftravelor one second's 
time for the car to reach equilibriumin rol] slip angles, and lateral ac- 
celeration The vaw accelerationappears as the difference between front 


andrear acceleration;it can be seento change sign during the transition 


indicating considerably less than critical yaw damping A second mode 
of oscillation witha frequency of approximately cps is also evident in 
tne slip angie and pat angle record 

These test data show conclusively that the directional control and 
stability of a caz ; very similar to that of a boat or an aircraft In all 
these cases, the directional control process appears as an oscillating 
system OI! at east tw legrees OI itreedadom that S,aSa proble m ol engi- 
neering mechar t is therefore of interest to see if equations can be 
develops d whic naes rive tnis actionana wt cl are Na S¢é don kr own proper 
ties ot tires susper Y teering ystem, and the car's inertia The 
Appendix shows thatsuch ananalysis in its simplest linearizedform leads 
to a fourtl yrder differential « juatior lr Ww ch one variable here the 
yaw velocity, is treatedas the dependentvariable, whilethe steering angle 
is used as theindependent variable Numerical values, corresponding to 
a representative medium siZ¢ ir, are used to compute irequency and 
damping of the two modes of free oscillatior These numerical values, 
however were not measured on the test car itself but are based on 
measurements On various similar cars 


The conditions of Case 1 in Table Icorrespondto the test run in Figure 


l and show fair agreement The lower mode has a calculated frequen- 
cy oO! 72 ps which is somewhat slower than that shown by the te st; its 
specific damping of .69 seems to fit the test quite well. The higher mode 
witha calculated frequency of 1.07 cpsand a specific damping of .38 does 
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kerboard Test with Step Input at Steering Wheel 














not correspond to the higher mode of about 2.0 cps indicated by the test 


This may be due to a wrong choice of the numerical values, but it is also 
possible thatthis apparent 2 cps motion is not real. It may simply repre- 
sent scatter of test data as stated in the original paps r> The two calcu- 
lated modesof .69 cpsand 1.07 cps caneasily combine so that inthis test 


they could not be separated 
Comparison of Case 1] and Case 2 shows that doubling the speed has 
only a minor effect on the two modes for this car, so that a driver would 


probably not detectany change This is usually true for well-behavedcars 


Both the test andthe analysis fall short of expressing directional sta- 


bility in a satisfactory way for two reasons They do not show the dis 
turbance produced by a road wave ora gust of wind This shortcoming 
could fairly easily be corrected by an extension of test and analysis 
Besides this, and more fundamental, is the objection that a driver's 
impression is based very much on steering wheel ''feel''; that means on 


the reaction forces he receives at the steering wheel during a maneuver 
and ona wavy or windy road We know that changes in caster and in the 


reversibility of the steering mechanism can greatly change whatthe driver 


calls directional stability The checkerboard test and the analysis of the 
Appendix are concerned with the properties of the ar under ''fixed con- 
trol'’ conditions, and hence contain no information onsteering reaction 
torces 

Some ofthe desiredinformation can be obtained by running a checker- 
board test under ''free control. ' This consists in throwing the car off 
r lan » Ay rking +} -+ ——t. sx) } ] na tt r r } ne T TI r ll 
Dalance DY j1erki! steering wheel andadtnen releasing it 1e Car Will 
then oscillate 1HDO0UtT a straight line while the steering whee swings back 

5 = 5S 

and fort} Such a maneuver is depicted in Figure 2 The frequency and 
decay ofthe car's oscillation can be observed without much instrumentation 
they correlate well with some items of steering performance in normal 


driving -for example with action on an S curve and with recovery from a 
x 


partial skid Such tests have been run on the checkerboard, but no com- 


pletely processea data are available 

The analysis canalso be modified to simulate roughly this typeof con- 
trol If the front wheels caster freely with zero friction in the steering 
linkage their cornering forces must be zero at all times; hence,setting 
the front whee yrnering power equal to zero, we get a first approxi- 
mation of the free control stability equatior 


ase 3 of the Appendix shows that the calculation tor this condition 


} Cs a 4 -l } . = as " . we : } 2A f } : 
leads to a first mode with a irequency ol 18 cps ana .5% speciiic Gamping. 
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, : : . . 
mon tne road usua shows trequencies nigner tnanthnis 


, Ssome- 


sh as 2 cps; however, the analysis correctly predicts a much 


ific damping for free control than for fixed control. Compari- 


of Case 3 and Case 4 shows that the Gamping ol both modes decreases 


Vitn ir reasing speed vwhichalso checks with road observation 
. } } ” te . ba ¢ . 

s will actually show a build-up of this oscillation at high car 

icating negative damping However, change of static friction in 


| | 5 | a Re 7 
ing linkage usually affects this free control oscillationmarkedly 


11a De 1 mistake todraw too many detailed conclusions from this 


with the addition of the free control test and free control analy- 


= an to 
neariver cCailis 


f the driver fully or partially prevents the wheel from moving 
bitrary maneuver. After we know this, we will of course ask 


he driver reactsto tne car s motions and torces; we will want to 


+ " | - + “ - ] * oc hi 
hat items inthis collection of data really correspond to his 
yn of directional stability and of responsiveness As a means of 


S tinal yrrelation Our past approacn 1s not very promising 


ill probably have to use techniques of servo motor analysis by 
ting a blockdiagram suchas Figure 3, according to the suggestion 


Chenea The input here is a disturbance produced by the 


f Vind, the ipper row I DIOCKS representing transter tunctions 
which establish the values of significant variables produced by 
tional properties of the car The driver is represented by the 
k which contains his transfer function or functions; he picks up 


im some variables, acts on them, and produces as output a new 


iT This is added to the original disturbance, thus forming a 
[ 5 rr 
esentwe do not know enoughabout the humanoperatoras a servo 
iown the transfer function we need inthis diagram; it is even 
that our present mathematics is incapable ii giving uS a satis- 
inctior But research now in process on human reaction and or 
ithematics may soon give usthe proper tools We shall then be 
nsider the stability ofthe whole system of car and driver so that 
yecify the directional properties of the car which will make i 
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Figure 2 - Free Control Action 
Schematic Diagram of Observations with Impulse Disturbance 
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APPENDIX 
SOLUTION OF EQUATIONS OF MOTION 


List of Symbols 


Course angle (yaw angle) 


y Path angle 
Steering angle 
20 Center slip angle (referred to roll axis) 
OC} Front slip angle 
a2 Rear slip angle 
8 Roll angle 
F Lateral inertia force at C.G. 
Fy Front cornering force 
Fo Rear cornering force 


use 


ying 


and 


gust 


SAE 


National Influenceson American Passenger Car Design, 


., Car Control Factorsand Their Measurements, 











Qu 


jr 








m = 


i * 


a = 


< 
u 


According to 


angle and 


may write 


F 


The testda 


this will be 


Mass of car 
Moment of inertia, around an axis parallel 
to the roll axis - through the C.G. 
Moment of inertia, around yaw axis 
Whee lbase 
Horizontal distance from C.G. to front wheels 
Horizontal distance from C.G. to rear wheels 
Distance of C.G. above roll axis 
Car Spe € d 
Lateral acceleration at C.G. 
Front cornering power 
Rear cornering power 
Front roll-steer ratio 
Rear roll-steer ratio 
Roll spring rate about the roll axis 
coefficient about the roll axis 
references | and 2 thereis a linear relation between slip 


ornering force for slip angles below six degrees, so that we 


Cx 


ta also show a variation of cornering power with vertical load; 


neglected in order to avoid non-linearity. To the slip angles 


produced by the tire must be added steering angles due to lateral load in 


the suspension 


nering pow 


Additional steerin 


be assume 


} 
] 


These will simply beincluded by modification of the cor- 
er © 


angles are produced by roll displacement and will 


d proportional to roll angle. This leads to 


uw So. ¢ 
l C) . 
F., 
=- +46 
ec 2 


or, solving for F: 


~ 
— 





_ 


From Figure 4 we can then write three equations of equilibrium. 


Lateral forces: 


Moments about yaw axis: 


Moments about roll axis: 
I 6+f£60+K®-mh)A=0 


Additional equations are derived from geometric relations based on 
Figure 5 
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™m 
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V, sind. - V. sin dG 
2 2 oO re) 
b 
For smallangles where dQ, sin Q, and Vv; z: V, 2: Vo = V, these three 


equations are: 


From (4b) and(4c) 


(1/L) |jads - b (y-4a)) Hence, 


“— 
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Figure 4 - Attitude of Car in Transient Turn 








5 - Relations Between Slip Angles at Front, Center 





, and Rear 
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(5) 


A last equation results from the consideration that the lateral acceler- 
ation at the C.G. is the sum of lateral acceleration at the roll axis plus 


lateral acceleration relative to the roll axis. 
L\=VW-hé (6) 


Considering V as constant, yas the independent variable, we have six 
dependent variables yw, 0, ©, A, dj, A>, so the six equations (1) to (6) 


are sufficient. We may solve them for one dependent variable, for ex- 


ample the yaw velocity 2, and find the equation 


~ 
~~ 


pv + a gpl + az 9+ ao + a) D = b3y+ bo Yt b, y+ boy ( 


Where the parameters have the values: 
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No complete setof 


numericalmeasurements for any one car is available, 


but values fora representative medium sizedcar, partly based on measure- 


ments and partly or 


With these valu: 
ft/sec (82 mph), the 
Cj = C2 14, 0Q0 


At OU itt é 
rT rel 
a, | 
4 
4 4 
as h 7 
. yORS 
ma 
For free contro 
At 60 ft/se 
i4 5o.5 
az 123 
ao 41% 
ay 2U93 


Hence the character 


estimates, would be as follows: 


000 lb/rad 
100 lb/rad 
500 lb ft/rad 


sec raqa 


speeds of V 60 ft/sec (41 mph) and V = 120 


ilues of the parameters are for fixed control (i.e., 


raqd) 


At 120 ft/sec 


AS a 


al Sita 

oie) 107.25 

69 340.22 
1163 


14,000 lb/rad 


At 120 ft/sec 


——— a a ke 


6.06 

3] 116.92 

35 248.61 
2093. 


{- . 
stic equations of equation ({/) are for 
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FIXED CONTROL: (i.e., C) = C2 = 14,000 lb/rad) 


Case (1) V = 60 ft/sec 


D* + 14.18 D> + 140.06 D“ + 676.69 D+ 2085 


Factors are: (D2 + 8.69D + 39.328)(D* + 5.49D + 53.024) 


= 0 


Roots are: -4.345 +4 4.522li and -2.745 + 6.74461 


Frequencies are: 2cps and 1.07 cps 
Specific dampings are: .69 and .38 


Case (2) V 120 ft/sec 


p* + 8.72 D2 + 107.25 D@ + 340.22 D+ 1163 =0 


Factors are: (D“ + 3.527D+ 15.936)(D2 + 5.193D + 72.999) : 


Roots are: -1.76354+3.5814i and -2.5965 + 8.1399i 


Frequencies are: .57 cps and 1.30 cps 
Specific dampings are: .44 and .30 
FREE CONTROL: (i.e., C, = 0, Co = 14,000 lb/rad) 
Case (3) V = 60 ft/sec 


D* + 8.87 D2 + 123.31 D@ + 418.35 D+ 2093 


y 


Factors are: (D* + 3.545D+ 27.042)(D2 + 5.325D + 77.390) = 


= 0 


Roots are: -1.7725 + 4.8888i and -2.6625 + 8.3846i 


Frequencies are: .78 cps and 1.33 cps 


Specific dampings are: .34 and .30 

















Case (4) V = 120 ft/sex 





D* + 6.06 D>? + 116.92 D“ + 248.61 D+ 2093 = 0 


Factors are (D¢ + 1|,.646D + 4 609)( D2 + 4.414D + 85.045) = Q 


= 0 Roots ar 823+ 4.8920i a 207 + 8.9540i 
f er € ire 8 Dp and 1.43 cps 
Sp amy are | ind 4 





Condition or lst Mode 2nd Mode 
if  ——_——————— me eee ————_—__——_—__F— 
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GENERATING MACHINES 


The principle of operationof a straight bevel gear generatoris shown 
inFigs. 3 to 5. There is a cradle carrying two reciprocating tools, each 
to generate one side of a tooth. A gear segment (Fig. 3) or a geartrain 
provides the constant "ratio of roll’ betweencradleand gear blank, which 
rotate aboutintersecting axes. The tool slides are set to a ''tooth angle’ 
to produce the required taper on tooth thickness (Fig. 5). 


The tooltips travelin a plane perpendicular to the cradle axis,which 
is the root plane of the bevel gear. The instant axis makes a smallangle 
with this plane. Thus, the ''generating gear'’ has less thana 90-degree 
pitch angle, which makes the profile curvature slightly greater than that 


generated by anexactcrown gear. Some newer straight bevel gear gener- 





ators allow adjustments for any pitch angle of the imaginary ''generating 


gear 


Figures 6 and 7 show theprinciple of aspiralbevel gear generator. 


CONJUGATE GEAR TEETH AND INTERCHANGEABLE GEARS 


Any tooth profile or lengthwise tooth formation of conjugate racks 
yields conjugate gears in the sense that for eachpoint of one tooth surface 
there is only one mating contact pointon the mating surface. Sang’ s! theo- 
ry, developed a hundredyearsago, emphasized the symmetry ofthe paths 
of contact on either side of the rack pitch line and for opposite sides of 
the teeth as a condition to obtain interchangeability. 


Fig. 8 shows the conjugate racks for interchangeable spur gears. Fig. 
9 shows conjugate gears in the form of an internal and a counterpart ex- 
ternal generating gear for conjugate spur gears; andalso shows the gener- 
ating gears of conjugate bevel gears. In this case, the pitch surfaces may 
»e planes (for crown-gear generation) or matching cones, and the tooth 





surfaces can be ofany convenient shape, unless interchangeability is re- 
quired. 


In many spiral bevel gears the generating surface of the cutter is a 
section of a cone. In the original conception, the axis of this cone (cutter 
axis in Fig. 10) is perpendicular to the pitch surface of the ''generating 
year. The ''generating gear'' lies ina plane tangent to the pitchconeat 
he instant axis. This type of cutter and its application was invented by 
James E. Gleason* some 40 years ago. In this case, the teeth are of 
constant depth. To obtain tapered tooth depth as in Fig. 6, certain de- 
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partures from exact theory are made, for which there is a practical so- 
lution. 


However, exactly conjugate teeth in bevel gears are not desirable 
under operating conditions, for the same reasons which have led more 


recently to profile modification and crowning of cylindricalinvolute gears. 


FORMATE* BEVEL GEARS 


A characteristic of any gear-generating method is that the shape of 
the generating tool 1s not reproducedon the surfaceof theteeth of the gear. 
Forming methods of gear cutting do reproduce the tool shape, and such 


methods are generally faster and can be performed on more rigidly built 
machines 


For many yearsit has beenmore economical tomake the larger mem- 


ber of a pair of spiral bevel or hypoidgears form-cut, rather than gener- 
ated. 





Then the mating pinion is generated using a cutter similar to the 
ee er 
form cut gear with which the pinionis to run. It is readily seen that a 


Formate gear and pinion can be made correctly conjugate in this way. 


SURFACE CHANGES 


Bevel gear generation permits more freedom than the processes in 
use for cylindricalgears. The tooth sides are often cut by separate tools 
or in separate operations, even in different machines, so that tooth thick- 
nessis not determined by tool dimensionand cutting depth. Modifications 
of pressure angle and spiral angle, even when one tool is employed, are 
used to control tooth contact or bearing. An interesting example is the 
use of a '' generating gear" with less than 90° pitch angle, to increase the 
profile curvature of the bevel gear by a few thousandths ofan inch to modify 
the contact with the mating gear. Spur gears could also be modified in 
this manner if similar tools and methods were in use. 


POINTS OF CONTACT, LINES OF CONTACT, AND CONJUGATE POINTS 


Eventhough conjugate bevel gears can be produced without exact knowl- 
edge of tooth profiles and lengthwise curvature, such knowledge is often 
desired to estimate tooth strength and surface durability. The methods 
usedare much the same as for spur gears and start withthe determination 


of points »f contact and lines of contact. 
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A line of contact between two meshing teeth is the line or curve in 





which two curved surfaces contact each other ata given instant. When 


motion occurs, this line changes in location, sometimes also in shape. 


A pointof such a line can be found by the fundamental principle (com- 
mon to spur and bevel gears but not tocrossedaxes gears) that the common 


normal ata point of contact passes through the instant axis. 


It is therefore only necessary to drop normals from the instant axis 
to the (instantaneously fixed) tooth surface in order to establish points of 
sucha line of contact. This procedure can be considerably simplified if 
the normals are dropped to the surface of the ''generating gear,'' whose 
shape is geometrically simpler than that of the operating gear, e.g., 
planein the crown gear of straight bevel gear teethand conicalin the crown 
gear of spiral bevel gear teeth. Lines of contact for some consecutive 


tooth positions determine the surface of action which, in spiral bevel gears, 





is a warped surface. Fortunately, we do not require these for general 
engineering purposes. 


The surface of the '' conjugate'' generated gear tooth can be established 





by transferring points of the projections of the consecutive common contact 


lines to a selected fixed position of the required conjugate gear tooth. 
J g } 


This drawing procedure is cumbersome because it is necessary to 
establish three different projections of each contact point and each conju- 
gate point. Locations of such points can be determined as accurately as 


desired by numerical calculations based on triangulation. 


TOOTH LAYOUTS 


The shape of a tooth section is necessary for estimating its strength. 
Cylindrical helical gears are not studied in a true section in a normal 
plane, but by substituting an equivalent spur gear as a satisfactorily close 





approximation. The equivalence is with respect to the size of the pitch 
circle, the circular pitch, circular tooth thickness, addendum and deden- 
dum, pressure angle, and profile shape. Similar approximations are 
satisfactory for bevel gears. 


TRANSVERSE SPHERE AND TRANSVERSE PLANE 


Profile curves taken on a spherical surface centered at the apex point 
will contact one another in the same manner as the transverse profile of 


spur gears. Figure 11 shows such transverse sphere with pitch circles, 
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profile curves, etc. Spherical geometry could be used to obtain exact 
results. Approximations based on plane line drawings are usually pre- 


ferred. 


Fora straight bevel gear, Figure 12, anequivalent spur gear is substi- 
tuted, drawn in a transverse plane perpendicular to a pitch cone element 
usually at the outer end of the tooth or, for engineering purposes, at the 
mean section. This method is known as Tredgold's> approximation and 
assumes the eguivalent spur gear on the surfaceof a back cone unwrapped 





into a plane. Its pitch radius is the back cone distance or the normal 
distance of a pointon the pitchcone elementfrom the gear axis. It is also 
4 


theorem for the curvatureof surfaces, it is the radius of curvature of the 


the radius ofa sphere inscribedinthe pitch cone. According to Meunier's 


pitch cone surface in the chosen transverse plane. 


NORMAL LAYOUT OF SPIRAL TEETH 


The pitch radius in the normal layout of a spiral bevel gear may be 
obtained by Euler's* theorem as the transverse radius divided by the square 


of the cosine of the spiral angle as graphically illustrated in Figure 13. 


The approximate tooth profiles of the bevel gear are obtained by roll- 
ing the equivalent pitch circle of the generating gear with that of the bevel 
gear. The generating crown gear is represented as a rack with straight 


pitch line. 


SUMMARY OF SPECIFICATIONS AND MACHINE SETTINGS 

While approximations are adequate for tooth layouts, exact calcu- 
lations must be made for cutter specifications and machine settings. Many 
different sizes and types of machines exist, and two or three generating 
methods may be used, even on the same machine. This does not allow 


standardization to the same degree as in spur and helical gears. 


INSPECTION AND MEASUREMENT 


The problem of bevel gear inspection is, in general, not approached 
by the same methods of exact measurement that are used for cylindrical 
involute gears, because the size and shape of the bevel gear tooth varies 
from toe to heel and becausethere is more than one way to setup a gener- 
ating machine for gears of given numbers of teeth and pitch. Also, there 


is no simple way to follow the tooth profile reliably and accurately. The 
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position at which a measurement is taken becomes as important as the 
measured dimensions; also the machine, the tools, the setup, and the 
selected generating motion affect the tooth thickness. There are more 
variables eventhan in worm gear teeth, which are notordinarily inspected 
for tooth shape. 


The real test of bevel gears consists of running a gearand pinion to- 
gether and observing sound and tooth bearing (Fig. 14). The so-called 
V-and-H check" is an example. In the central relative position the gears 
should contact withonly a portion of the tooth surface. The tooth bearing 
canbe shifted from one end of a tooth to the other by slight changes in the 
relative position of the gears; and the required amount of change serves 


as a measurement of the ''ease-off.'' One setof bevel gears may be used 
as masters or ''gages'' to check other production gears. This running 


test quickly and easily assures quality and duplication for interchange- 


ability, which is the ultimate goal of gear manufacture. 
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